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Abstract 

Hot interstellar gas in elliptical galaxies has two 
sources: mass lost from evolving stars and a much 
older component that accompanied galaxy formation 
or arrived subsequently by secondary cosmic infall to¬ 
ward the galaxy group containing the elliptical. We 
present here an approximate but comprehensive study 
of the dynamical evolution of the hot gas in mas¬ 
sive elliptical galaxies born into a simple flat uni¬ 
verse. Baryonic and dark matter are both conserved. 
We use NGC 4472 as a prototypical massive ellipti¬ 
cal having a well-observed hot interstellar medium. 
We allow for star formation in a simple single burst 
using a Salpeter IMF but treat the gas dynamics in 
detail. The galaxy has a de Vaucouleurs stellar core 
and a Navarro-Frenk-White dark halo surrounded by 
inflowing cosmic matter. 

Using rather standard assumptions and parame¬ 
ters, we are able to successfully reproduce the gas 
density and temperature distributions - n(r) and T(r) 
- in the hot interstellar gas determined by recent X- 
ray observations. Our model is sensitive to the baryon 
fraction of the universe, the Type II supernova energy 
released per unit stellar mass, and the time of galaxy 
formation. But there is some degeneracy; as each of 
these parameters is varied, the effect on model fits 
to n(r) and T(r) is similar. Nevertheless, secondary 
inflow of cosmic gas is essential for successful fits to 
n(r) and T(r). 

Some gas is expelled from the stellar galactic core 
at early times when the Type II supernova energy was 
released. As a result, the present day baryonic frac¬ 
tion has a deep minimum in the outer galactic halo. 


Interstellar gas that cooled since the time of maxi¬ 
mum star formation cannot have all collected at the 
galactic center but must be spatially dispersed; oth¬ 
erwise both gas temperatures and stellar dispersions 
in the galactic center would be larger than those ob¬ 
served. 

Finally, when relatively gas-rich, X-ray luminous 
models are spatially truncated at early times, sim¬ 
ulating tidal events that may have occurred during 
galaxy group dynamics, the current locus of truncated 
models lies just along the L x , X-ray size correlation 
among well-observed ellipticals. This is another strik¬ 
ing confirmation of our model of elliptical evolution. 
Subject headings: galaxies: elliptical and lenticular - 
galaxies: formation - galaxies: evolution - galaxies: cool¬ 
ing flows - X-rays: galaxies 

1. INTRODUCTION 

In this paper we present a new, more comprehen¬ 
sive model for the evolution of hot gas in massive 
ellipticals and dominant ellipticals in galaxy groups. 
For the first time we combine gas produced by stel¬ 
lar mass loss with additional gas that flows into the 
galaxy regarded as a perturbation in the Hubble flow. 
The radial distribution of gas density and tempera¬ 
ture in bright ellipticals determined by recent X-ray 
observations can be understood only by combining 
these two sources of gas, stellar and cosmic. 

Traditionally, it has been assumed that most or all 
of the hot gas observed in ellipticals arises as a result 
of mass loss from the aging stellar population, but re¬ 
cent observations with the ROSAT satellite indicate 
a more complicated origin for the gas. Davis & White 
(1996) showed that the temperature of the hot gas is 
generally ~ 1.5 times hotter than the equivalent tem¬ 
perature of the stars T*. This implies that the gas is 
in virial equilibrium in a deeper, and therefore larger, 
potential than the stars. Even more relevant, X-ray 
images from ROSAT have revealed very extended re¬ 
gions of X-ray emission surrounding many massive 
ellipticals, extending far beyond the optical image. 
Mathews and Brighenti (1998) have shown that the 
emission from these extended regions can dominate 
the overall X-ray luminosity. This realization resolves 
a long-standing puzzle for elliptical galaxies: the enor¬ 
mous variation in X-ray luminosity L x in ellipticals 
having similar optical luminosity. However, the total 
mass of hot gas in the outer regions of many of these 
galaxies is too large to be understood in terms of stel- 
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lar mass loss alone. In Brighenti & Mathews (1998) 
we showed that agreement with the observed hot gas 
density and temperature profiles is greatly improved 
by assuming that a large mass of “circumgalactic” gas 
existed around the galaxy at early times. Much of the 
gas that filled the galaxy at the time of galaxy forma¬ 
tion is still present today. This raises the possibility 
that the hot gas in elliptical galaxies can reveal, like 
growth rings in a tree, the sequence of events that oc¬ 
curred long ago during the epoch of galaxy formation. 

In this paper we extend this idea to consider the 
continuous inflow of gas into the galactic halo over 
the Hubble time. This type of flow is a natural con¬ 
sequence of the large perturbation in the Hubble flow 
that formed the galaxy and the small group in which 
it may reside. Our approach in constructing these 
models is to adopt the simplest or most plausible 
procedure. We assume for example that the evolu¬ 
tion of the dark matter can be described by the nor¬ 
mal self-similar flow in a perturbed flat universe, but 
as the central stationary dark halo grows in size, it 
adapts to a shape dictated by more detailed N-body 
calculations. The baryonic gaseous component flows 
in with the dark matter but eventually shocks and 
compresses to the virial temperature ~ 10 7 K of the 
dark halo. From the very beginning, as X-ray en¬ 
ergy is radiated away, cold gas condenses at the cen¬ 
ter of the flow. When a sufficient mass of baryons 
has accumulated, we form the stellar distribution of 
the elliptical in a single event. In this way we circum¬ 
vent the complex dynamical and merging processes 
that occurred during the formation of the big ellipti¬ 
cal. At the same instant when the galaxy forms the 
remaining gas near the galaxy is heated by Type II 
supernovae and the galactic stars continue to evolve 
as a single burst thereafter. In spite of the simplic¬ 
ity of this model, the subsequent evolution of the gas, 
both that expelled from the stars and gas continu¬ 
ously flowing in toward the perturbation, is treated 
in a self-consistent, reasonably realistic manner. This 
is confirmed by the agreement of our results with cur¬ 
rent X-ray observations. It is gratifying that we can 
match so well the observed mass and entropy distri¬ 
bution in large ellipticals with our models and still 
conform to global cosmological and astrophysical con¬ 
straints set by the conservation of baryonic mass, dark 
mass and supernova energy. 

In this study we are concerned with establishing 
the global properties of hot gas in ellipticals as a con¬ 
sequence of gas dynamical evolution from the earliest 


times. Therefore, in order to avoid unnecessary de¬ 
tail, we shall not discuss here the detailed metal abun¬ 
dance in the gas nor shall we consider the important 
effects of galactic rotation. 

Although our results are generally applicable to all 
massive ellipticals and the group environment they in¬ 
habit, as a rigorous test of the general success of our 
models it is useful to use parameters specific to a par¬ 
ticular well-observed galaxy; for this purpose we have 
chosen NGC 4472, the central dominant elliptical in 
a subcluster in the Virgo cluster. 

2. PROPERTIES OF NGC 4472, A LUMI¬ 
NOUS ELLIPTICAL 

The X-ray surface brightness variation in massive 
ellipticals can be inverted to determine the spatial 
variation of electron density with physical radius. The 
radial plasma temperature profile can be similarly de¬ 
termined from variations in the apparent X-ray spec¬ 
trum viewed in projection. Since the hot interstel¬ 
lar gas is in near hydrostatic equilibrium, the total 
gravitating mass M(r) can be found from the tem¬ 
perature and density profiles. At present the best 
galaxies for this purpose are those relatively nearby 
and luminous ellipticals that have been imaged by 
both Einstein HRI and ROSAT HRI & PSPC: NGC 
4472, 4649, and 4636 (Brighenti & Mathews 1997b). 
Of these we chose NGC 4472 for comparison with our 
gas dynamical models since the X-ray data for NGC 
4472 exists over the largest range, ~ 0.016 to ~ 16 
effective radii where r e = 8.57 kpc at a distance of 
17 Mpc. Brighenti & Mathews (1997b) showed that 
the hot gas in NGC 4472 is at rest in the stellar po¬ 
tential in the range 0.1 — 1 r e where the total mass 
(assuming hydrostatic equilibrium) is identical to the 
stellar mass predicted from the dynamical mass to 
light ratio based on a two-integral stellar distribution 
function (van der Marel 1991). This excellent agree¬ 
ment indicates that gas pressure alone is sufficient to 
support the hot gas; no additional pressure (magnetic, 
cosmic rays, turbulence, etc.) is needed nor indicated 
in 0.1 — 1 r e . Unfortunately, beyond a radius of ~3.5 
r e the X-ray image of NGC 4472 becomes asymmetric 
probably as a result of an environmental ram pressure 
interaction or an interaction with the dwarf irregu¬ 
lar galaxy UGC 7636 (Irwin & Sarazin 1996; Irwin, 
Frayer, & Sarazin 1997). However, the azimuthally 
averaged X-ray surface brightness used by these au¬ 
thors to determine the average density profile is glob- 
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ally representative of luminous ellipticals. NGC 4472 
is a prototypical bright elliptical that exhibits all of 
the characteristics of these massive galaxies. 

The upper panel of Figure 1 shows the gas den¬ 
sity (Trinchieri, Fabbiano & Canizares 1986; Irwin 
& Sarazin 1996) and temperature (Irwin and Sarazin 
1996) profiles in NGC 4472 together with analytic fits 
as described by Brighenti & Mathews (1997b). The 
corresponding total mass M (r) is shown in the lower 
panel of Figure 1 together with best fitting stellar 
and dark halo mass distributions. The stellar mass 
is described with a de Vaucouleurs profile with effec¬ 
tive radius r e = 8.57 kpc. The stellar luminosity is 
L b = 7.89 x 10 10 L bq corresponding to a total stel¬ 
lar mass M* = 7.26 x 10 11 M 0 using a mass to light 
ratio of AI*/L B = 9.20 (van cler Marel 1991). The 
dark halo in NGC 4472 dominates the total mass for 
r > r e . As shown in Figure 1 a dark halo described 
by the Navarro, Frenk & White (1996) density profile 
(NFW halo) using a total halo mass Mh = 4 x 10 13 
M 0 fits the dark matter distribution quite well at 
large radii. [A reservation about the applicability of 
the NFW profile to NGC 4472 expressed by Brighenti 
& Mathews (1997b) is now seen to be in error.] The 
density distribution in an NFW halo is described by 

_ 200 c 3 p H 

Ph 3 ln(l + c) — c/(l + c) (r/r s )[l + (r/r s )] 2 

where ph is the current mean density of the universe, 
r s = Tg n /c and r gn is the current virial radius of the 
halo where the density is ~ 200 p B - The concentration 
parameter, c = 10.72(M/,/3.3 x 1O 13 M 0 ) -0 1233 , is 
based on Figure 8 of NFW. All galactic parameters 
are scaled to our adopted distance to NGC 4472, D = 
17 Mpc. 

The agreement of the mass distribution of NGC 
4472 with an NFW halo plus a de Vaucouleurs core 
is good but not perfect, particularly at the transi¬ 
tion region near ~ 10 kpc. The shape of the inner 
part of galaxy-sized dark halos is controversial, e.g. 
Kravtsov et al. (1998) predict a flatter density profile 
within the NFW characteristic radius r s ; such a halo 
would improve the fit to the total mass in Figure 1. 
However, any mass profile determined by N-body cal¬ 
culations in the absence of baryons will be centrally 
steepened and redistributed somewhat by the dissi¬ 
pative inward migration of the baryon component. 
Nevertheless, in the calculations described here we 
shall adopt the combined de Vaucouleurs and NFW 
mass distribution for NGC 4472. The relatively small 


discrepancy apparent in Figure 1 has essentially no 
significant influence on our gas dynamical models. 

We have taken published profiles of the electron 
density n(r) directly from the literature although 
Trinchieri, Fabbiano & Canizares (1986) and Irwin 
& Sarazin (1996) have made different assumptions 
about T{r) and the metallicity variation z{r) in the 
gas. Since the indicated gas density depends on T(r) 
and z(r) as they influence the cooling rate, n(r) oc 
A(T, 2 ) -1 / 2 , the “observed” values of the gas density 
in Figure 1 could change somewhat when more defini¬ 
tive information on T(r) and z(r) becomes available. 


3. EVOLUTION OF THE DARK HALO BY 
SECONDARY INFALL 

For simplicity in estimating the growth of the dark 
halo in NGC 4472, we assume a critical = 1 uni¬ 
verse and a current age t n = 13 Gyrs corresponding 
to H = 50 km s -1 Mpc -1 . The self-similar flow of 
dark matter toward an overdensity perturbation in a 
critical universe (Bertschinger 1985) consists of two 
regions: a smooth outer flow in which the Hubble 
flow reverses at a turnaround radius rt a (t) and an in¬ 
ner flow consisting of an ensemble of caustics or den¬ 
sity cusps produced by orbital wrapping of the col¬ 
lisionless dark fluid as it oscillates about the center 
of the perturbation. As infalling dark matter accu¬ 
mulates in the dark halo, its mass is added to the 
outer part, leaving the central mass distribution es¬ 
sentially undisturbed. Although newly arrived shells 
of dark matter pass through the central regions of 
the previously existing dark halo, they move very 
rapidly there and do not contribute appreciably to 
the time-averaged density except near the outermost 
caustic. Averaging over the caustic peaks, the total 
dark mass density in Bertschinger’s similarity solution 
varies as pd oc r -9 / 4 , but this is much too centrally 
peaked when compared to typical N-body simulations 
(Navarro, Frenk & White 1996). 

We resolve this difficulty by assuming that an 
NFW halo forms as a stationary core at the center 
of an outer smooth flow described by Bertschinger’s 
similarity solution. The outer flow must be attached 
to the stationary NFW halo so as to conserve total 
dark mass at all times. In terms of the similarity 
variable 

r 

r ta (t) 
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the outer solution can be fit quite precisely with 


M( A) = A 3 


a + b( A/Ai) 2 
1 + (A/Ai) 2 


where a = 3.4970, b = 5.8981 and Ai = 1.1296. This 
analytic fit to the smooth outer flow onto a point mass 
is exact at A = 0, 1 and oo. 

The current virial radius r gn within the self-similar 
dark halo corresponds to an overdensity of ~ 200, i.e. 

( 3 M n V /3 

rgn ^ 4^200 p Hn ) 


where pn n — ^-I^Gt 2 n is the current total cosmic 
density and M n = 4 x 10 13 Mq is taken to be 
the current virial mass. In the self-similar solution 
the approximately stationary halo confronts the sec¬ 
ondary inflow at the radius of the outermost caustic 
r cn = 1.528r 9 „ at time t n . The current turnaround 
radius rtan = Tt a {t n ) can be found by insisting that 
the outer self-similar solution meet the NFW halo at 
radius r cn where the interior masses of the two mass 
distributions must agree: 

4 .3 A At / \ _ n r f( r cn/ r s) 

Q ^ tanPHn-MyTcn/Ttan) Mn r/ \ 

3 /(c) 

where 

/( x) = ln(l + x) - — 

1 + x 

describes the NFW mass profile. With our cho¬ 
sen parameters, the current turnaround radius is 
rtan = 1.096 x 10 25 cm. By similarity arguments, 
the turnaround radius at earlier times is rt a (t) = 

Ttan(t/t n ) 8/9 . 

Finally, the radius r c at which the NFW halo meets 
the secondary inflow is found by solving 

^[r ta (t)} 3 p H (t)M[r c (t)/r ta {t)} = M n ^ 

for r c (t), keeping r s constant. Although A4(A) is 
unvarying with time in the similarity solution, we 
have relaxed this condition in the expression above 
to achieve a fit to the NFW halo. The total mass dis¬ 
tribution in the stationary halo and the external flow 
are illustrated at several times in Figure 2. 


4. BARYONS AND SINGLE BURST STAR 
FORMATION 

Within the context of the 0 = 1 cosmology the 
baryonic contribution to the total density 0& is con¬ 
strained by primordial nucleosynthesis to a narrow 


range, 

0.04 ^ 0 & (R/50) 2 £ 0.06 

(Walker et al. 1991). For a given value of O^ the mass 
distribution in Figure 2, multiplied by (1 —0&), repre¬ 
sents the varying potential of dark matter. Provided 
the temperature of the baryonic gas is much less than 
the virial temperature of the dark halo, T < 10 6 K, 
cosmic gas will flow inward toward the galaxy pre¬ 
cisely following the dark matter until the gas encoun¬ 
ters a shock front. Within this accretion shock the 
mass distribution of the baryonic component is deter¬ 
mined by gas dynamics and the assumed model for 
star formation. 

Interior to the outward facing shock the cosmic gas 
is heated to T ~ 10 7 K, similar to the virial tempera¬ 
ture of the dark halo. After a few Gyrs a large mass of 
gas will have collected and radiatively cooled within 
the dark matter halo; it is natural to assume that this 
gas forms the old stellar system observed today. In 
our very simple model for star formation we assume 
that all the stars in NGC 4472 form as a single burst 
at time i*. One possible choice of f* would be to as¬ 
sume that stars form from the cooled gas as as soon 
as a mass of cold gas equals the current stellar mass of 
NGC 4472. However, dynamic considerations suggest 
that the de Vaucouleurs stellar profile results from 
the merging of stellar systems that must have formed 
prior to the merging events, so stars are likely to have 
predated NGC 4472. Furthermore, gas cools more ef¬ 
ficiently and can form stars somewhat sooner in an 
earlier generation of small galaxies in which the dark 
matter and cooling gas are more concentrated than in 
NGC 4472. For these reasons we regard the time of 
star formation f* as a variable parameter in our mod¬ 
els; if the total mass of radiatively cooled gas is less 
than the mass M* of NGC 4472 at time t*, some of 
the hot gas is used to make up the deficiency since it 
could have cooled by this time in smaller pre-galactic 
stellar systems. At the time of star formation the gas 
density is reduced by a factor [M g (r s h) — M*\/M g {r s h) 
where M g (r s h) is the total (all of the cold and per¬ 
haps some hot) gas mass within the shock radius r s h 
just before f*. Experience has shown that this arti¬ 
ficial mass rearrangement at f* results in a relatively 
short-lived gas dynamical transient. 

Soon after time t* massive stars (m > 8 Mq) begin 
to produce SNII explosions which heat the remaining 
gas considerably. We deposit this SNII energy instan¬ 
taneously at time f* by increasing the specific thermal 
energy density e(r) by an appropriate amount every- 
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where within the shock radius at t*. The possibility 
that some SNII events predated the physical forma¬ 
tion of NGC 4472 will not be considered here. The 
combined energy of the SNII explosions is sufficiently 
large to drive an outflow from the galaxy and the 
associated galaxy group, creating a second outward 
moving shock that moves out beyond the accretion 
shock described above. 

As soon as stars form they begin to lose mass by 
normal stellar evolution. Most of this stellar mass loss 
occurs at early times shortly after f*. Unlike more re¬ 
alistic, rotating galaxies in which radiatively cooled 
gas forms into disks of size Tdisk r e (Brighenti & 
Mathews 1996; 1997a), in the purely spherical mod¬ 
els discussed here cold gas can only form at the very 
center. From f* until the present time enough cold 
gas can form at the center of the galaxy to influence 
the gravitational potential there, causing the gas tem¬ 
perature within r e to rise considerably above values 
observed. We avoid this possibility here by assum¬ 
ing that the total mass of gas that cools after time f* 
also forms into stars distributed at all radii within the 
stellar galaxy. The mass of stars that is assumed to 
form in this diffuse manner after time f* is much less 
than the current stellar mass of NGC 4472 and can 
be ignored in evaluating the global galactic potential. 
Therefore to implement this assumption, we simply 
ignore the self gravity of the cooled gas at the galactic 
center. We also ignore the relatively small SNII con¬ 
tribution resulting from this type of additional star 
formation. 

Although our model for star formation is obviously 
ad hoc and oversimplified, we believe that the subse¬ 
quent dynamics of the hot gas, which is our main 
interest here, is modeled reasonably well. 

5. GAS DYNAMICS 

The standard gas dynamics equations that de¬ 
scribe the evolution of hot interstellar gas in ellipticals 
are the usual conservation equations with additional 
source and sink terms: 


where e = 3kT/2pm p is the specific thermal energy; 
p = 0.62 is the molecular weight and m p the proton 
mass. M to t(r) is the total mass of stars, dark matter 
and hot gas within radius r. 

All three equations involve source terms describing 
the rate op* (gm s _1 cm -3 ) that gas is ejected from 
evolving stars and supernovae, i.e. a = a* + a sn . To 
estimate cc* (f) we assume a single burst of star forma¬ 
tion with a power law initial mass function (IMF) be¬ 
tween lower and upper mass limits mp and m u . Using 
the procedure described by Mathews (1989) we find 

a*( t) = a* (t n - f*) [t/(t n - i*)] -s(t) 

where t n = 13 Gyrs is the current time. For a Salpeter 
IMF of slope x = 1.35 with mp = 0.08 and m u = 100 
Mq we find a*(f„ — f*) = 5.6 x 10 _2 ° s _1 and s(t) 
varies slowly with time (s = 0.83, 0.97, 1.11, 1.34 at 
0.001, 0.1, 1, 11 Gyrs respectively). 

In the thermal energy equation A(T, z) is the op¬ 
tically thin radiative cooling coefficient which is a 
function of temperature and metallicity (Sutherland 
& Dopita 1993). Thermal energy lost from the gas by 
radiation is largely restored by the gravitational com¬ 
pression of the gas as it flows deeper into the galaxy. 
Radiative cooling eventually dominates in the densest 
parts of the flow. In the computed models we allow 
the gas to cool to 10 K. In addition we assume that the 
interstellar gas is heated by the dissipation of the or¬ 
bital energy of mass-losing stars and by supernovae. 
The mean gas injection energy is e 0 = 3fcT 0 /2p,m p 
where T a = (a*T* + a sn T sn )/a. The stellar tempera¬ 
ture T*(r) is found by solving the Jeans equation for 
a de Vaucouleurs profile with additional dark mat¬ 
ter beyond r e . Type la supernova heating is assumed 
to be distributed smoothly in the gas; we ignore the 
buoyancy of gas locally heated by supernovae that 
may transport energy and metals to larger galactic 
radii. The heating by supernovae is described by mul¬ 
tiplying the specific mass loss rate of supernovae a sn 
by the characteristic temperature of the ejected mass 
m S m T sn = m p E sn /2>km sn , i.e. 


dp 

dt 


l_ d_ 

r 2 dr 


(r 2 pu) = ap *, 


du 

~dt 


du 

l dr 


and 


dP GM tot (r) 

o P o ap*u, 

or 



P dp 
p dt 


p 2 A 

- 5 —b 0/3* 

m 2 


So 


— £ — 


P 

P 



a sn T sn = 2.13 x 1CT 8 SNu(i) (U sn /10 51 ergs) 

hr 1 - 7 (l b /l B q )~ 0 35 k s- 1 

where h = H /100 = 0.5 is the reduced Hubble con¬ 
stant and we adopt E sn = 10 51 ergs as the typical en¬ 
ergy released in both Type la and Type II supernovae. 
The Type la supernova rate is represented in SNu 
units, the number of supernovae in 100 years expected 
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from stars of total luminosity 10 10 Lbq. The depen¬ 
dence a sn oc Lg 0 35 derives from the mass to light 
ratio for ellipticals given by van der Marel (1991), 
M^t/L b = 2.98 x 10 ~ 3 L^ 35 h 17 in solar units. 

The equations are solved in one dimensional spher¬ 
ical symmetry using an appropriately modified ver¬ 
sion of the Eulerian code ZEUS. Spatial zones in¬ 
crease logarithmically and extend far out beyond the 
turnaround radius where the cosmic gas is superson¬ 
ically receding from the galaxy. Therefore we can 
assume “outflow” boundary conditions at the outer¬ 
most spherical zone, confident that no disturbance at 
the boundary can propagate back into the galactic 
flow. For our purposes here it is not necessary to con¬ 
sider the detailed thermal physics of the intergalactic 
gas; we simply assume that the temperature of this 
gas is isothermal at T = 10 K until it passes through 
the accretion shock. Although the intergalactic tem¬ 
perature is expected to be greater than this (due to 
photoelectric heating by quasar radiation for exam¬ 
ple), the value of the intergalactic gas temperature 
makes no difference in our calculations as long as it is 
small compared to the virial temperature of the dark 
halo of NGC 4472, i.e. < 10 6 K. 

6. SUPERNOVA HEATING 

Supernovae in ellipticals today are infrequent and 
are all of Type la. Cappellaro et al. (1997) find a cur¬ 
rent rate of 0.066 ± 0.027(i7/50) 2 in SNu units. The 
frequency of SNIa in the past is currently unknown 
although it seems reasonable to assume that it has 
been decreasing over cosmic time, 

SNu(t) = SNu(f„)(t/;f n )- p . 

Ciotti et al. (1991) recognized that the relative rates 
of stellar mass ejection (~ t ~ s ) and Type la super¬ 
nova (~ t~ p ) determine the gas dynamical history of 
hot interstellar gas in ellipticals. For example if p > s 
then galactic winds may occur at early times but if 
p < s winds may occur at late times. However, in 
the presence of large amounts of additional cosmic 
gas flowing into the galaxy, galactic winds driven by 
Type la supernovae will tend to be suppressed. We 
assume here that p = 1 and SNu (t n ) = 0.03 so Type 
la supernovae are unable to generate winds by them¬ 
selves. 

However, the enormous energy released by super¬ 
novae of Type II produced by massive stars at very 
early times can drive an outflow from the galaxy even 


in the presence of secondary infall. With a single 
burst, power law IMF described by 

4>{m)dm = 4> 0 m^ < - 1+x ^dm nig < m < m u 
the number of SNII per unit stellar mass formed is 

Nu x — 1 m$ x — m~ x 
1111 Mm x m]~ x - m\~ x 

where M* t is the total initial mass of stars. All stars 
greater than eight solar masses (m > m%) are as¬ 
sumed to become SNII. For the Salpeter IMF dis¬ 
cussed above (x = 1.35, me = 0.08, m u = 100) we 
find r/ji = 6.81 x 10 -3 SNII per M 0 . The specific su¬ 
pernova number r]u is sensitive both to the IMF slope 
and the mass limits; for example pjj = 11.6 x 10 -3 
for x,me,m u = 1.35,0.3,100 and pn = 16.1 x 10 -3 
for x,me,m u = 1., 0.08,100. We conclude that pn 
is uncertain by a factor of at least 2 or 3. If pn = 
6.81 x 10 -3 , the total amount of SNII energy pro¬ 
duced during the single burst formation of NGC 4472 
{M^t = 7-26 x 10 11 Mq) is 4.9 x 10 60 ergs assuming 
E sn = 10 51 ergs per supernova. An additional un¬ 
certainty is the efficiency e sn of communicating this 
SNII energy to the hot interstellar gas; if massive stars 
form near cold, dusty gas clouds a considerably frac¬ 
tion of SNII energy could be radiated away by dust 
and dense gas. In galaxy formation simulations it is 
assumed that e sn ~ 0.1 — 0.2 of the supernova energy 
is returned to the interstellar gas (e.g. Kauffmann, 
White & Guiderdoni 1993), but in general e sn must 
be regarded as an adjustable parameter. 

7. THE STANDARD MODEL 

Our objective is to solve the gas dynamical equa¬ 
tions from the earliest times to the present (t n = 13 
Gyrs) and identify model parameters that result in 
gas density and temperature profiles similar to those 
observed today. While the number of uncertain pa¬ 
rameters is formidable - H, fi*, t n , IMF, A (T,z), 
SNII energy, e sn , etc. - there are also robust con¬ 
straints set by cosmological density variations and the 
global conservation of dark mass, baryonic mass and 
energy. 

We begin by describing the results of our “stan¬ 
dard” model. The standard model is not chosen for 
its perfect agreement with n(r) and T(r) observed in 
NGC 4472, although that agreement is quite satis¬ 
factory, but because the parameters and procedures 
used in the standard model are the most plausible 
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and least controversial. In our standard model we as¬ 
sume f If, = 0.05 and form the galaxy at time f* = 2 
Gyrs. A Salpeter IMF is used with me = 0.08 and 
m u = 100 for which a*(t s ) = 4.7 x 10 -20 (f/f s ) -1 ’ 26 
s^ 1 is an approximate, time-averaged value of the stel¬ 
lar mass loss rate and t s = t n — f* is the current age 
of the stars. The IMF also determines the specific en¬ 
ergy generated by SNII explosions, rj st d = 6.81 x 10 -3 
SNII per Mq, and we assume that this energy is fully 
converted to the hot gas, e sn = 1. Clearly the impor¬ 
tant parameter for any model is r/iie sn . The current 
rate of Type la supernovae is SNu(t n ) = 0.03 and it 
has decreased with time as t~ p with p = 1. In the 
standard model the radiative cooling was varied with 
radius to simulate an increasing gas metallicity closer 
to the center of the galaxy, A(T, r) = A(T)/(r) where 
/(r) = max[0.3, 2 —1.7(r/r 2 )] with r 2 = 40 kpc. Such 
a variation is consistent with the expected radial vari¬ 
ation of gaseous metallicity and the metal-dependent 
cooling coefficients described by Sutherland & Dopita 
(1993). The function /(r) is designed to match a ra¬ 
dial abundance variation similar to that observed in 
the hot interstellar gas in NGC 4472 (Forman et al. 
1993; Matsushita 1997). 

In Figure 3 we compare the computed gas den¬ 
sity n(r, t n ) and temperature T(r, t n ) for the standard 
model with observations of NGC 4472. The fit to the 
observed temperature is excellent and the the density 
agreement is acceptable, within a factor of 2 of the 
observed values over a range of 100 in radius. Recall 
that the “observed” density variation is also some¬ 
what uncertain as discussed earlier. Nevertheless, it 
appears that the gas density observed beyond 10 kpc 
is higher than that of the standard model while the 
reverse is true within 10 kpc. We have found that 
the first of these problems can be corrected while the 
second is more difficult. 

Figure 4 shows in more detail the spatial varia¬ 
tion of gas density, temperature, pressure and velocity 
over the full range of our calculations at three times 
including t n . It is remarkable how little the central 
flow profiles have changed since t = 4 Gyrs. The frac¬ 
tion of the gas density within 100 kpc due to stellar 
mass loss decreases slowly with a*(<). While new gas 
continues to flow into the galaxy at later times, this 
gas accumulates largely in the outer regions which 
grow in size while the solution within ~ 100 kpc re¬ 
mains in a quasi-steady state. Notice that there are 
two shocks. At time t n the main accretion shock is at 
log Vkpc = 3 and it is preceded by a second, smaller 


shock at log?’fc pc = 3.4. Beyond this outer shock, at 
the turnaround radius, log?y an = 3.6, the perturbed 
Hubble flow velocity vanishes. The influence of the 
outer shock can be clearly seen in all flow variables. 
This outer shock begins at time t* when the SNII 
energy is released; as it moves outward it creates a 
transient baryonic outflow that pushes some gas far 
out into the galactic halo and beyond. At later times 
radiative losses just behind the outer shock are signif¬ 
icant. The small increase in gas density before enter¬ 
ing this outer shock is an expected feature of the per¬ 
turbed Hubble flow. This calculation was performed 
with an intergalactic gas temperature of 10 K, but the 
flow pattern within the outer shock was unchanged 
when the calculation was repeated with intergalactic 
temperature increased to 10 4 K. 

The radial variation of the local baryonic fraction 
at time t n is shown in Figure 5 for the standard solu¬ 
tion; the hot gaseous contribution to the baryon frac¬ 
tion is shown by the dashed line. Interior to ~ 100 
kpc most of the baryons are stellar, but the gaseous 
contribution becomes significant at 1 Mpc and domi¬ 
nant beyond 10 Mpc. The most interesting feature in 
this plot is the large baryonic minimum centered at 
~ 500 kpc - this feature is a relic of the gaseous evac¬ 
uation created by the release of SNII energy. Baryons 
formerly in this region were pushed outward, forming 
a region of baryonic overdensity (f lb > 0.05) just be¬ 
yond a radius of 2.5 Mpc. Clearly such a large radius 
is irrelevant for NGC 4472 since it would enclose the 
entire Virgo cluster. However, the baryonic depletion 
is real and the rapid spatial variation of f2&/$7 in the 
outer galactic halo must be considered in interpreting 
observed values of the baryon fraction. For example, 
ROSAT observations in NGC 4472 extend only to ra¬ 
dius ~ 150 kpc where the local baryon fraction is seen 
to be very steeply declining in Figure 5. Clearly, mea¬ 
surements of the baryon fraction at such galactic radii 
could easily be greater or less than the true cosmic 
value (Clb/Cl = 0.05 for our standard model). 

We show in Figure 6 the time variation of the mass 
of hot gas, cooled gas and stars in the standard model. 
At time f* = 2 Gyrs all of the cold gas and some of 
the hot gas is consumed to form the single burst of 
stars in NGC 4472. Thereafter the stellar mass is 
assumed to be constant, M*(t*) = M*(t n ); we have 
ignored the small (~ 20 percent) decrease in M* ex¬ 
pected from stellar mass loss since time t n . Although 
hot gas (T > 10 6 K) continues to increase within the 
main accretion shock, the amount of hot gas within 
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150 kpc is constant for t > 3 Gyrs, suggesting a 
quasi-steady state in this part of the flow. Cold gas 
continues to form after time f* reaching a final mass 
Mcg(t n ) = 4.8 x 10 10 Mq at the end of the calcula¬ 
tion. While M cg (t n ) is comparable to the total mass 
of hot gas within 150 kpc, it is only ~ 7 percent of 
the current stellar mass. Unlike our earlier models in 
which we assumed that all of the hot interstellar gas is 
generated by stellar mass loss, the additional supply 
of cosmic gas at later times in the standard solution 
also results in more cold gas overall. 

In computing the gravitational potential we have 
ignored the mass of cold gas, assuming that some 
or all of it contributed to later generations of stars, 
formed throughout a significant volume of the galaxy. 
However, it is interesting to relax this assumption and 
include the gravity of cold gas which, in our idealized 
non-rotating models, forms a concentrated point mass 
at the very center of the galaxy. This variant of the 
standard model is shown in Figure 7. Compared to 
Figure 3, the density profile in Figure 7 is slightly 
improved within a few kpc, but gas within about 5 
kpc of the galactic core in Figure 7 becomes heated 
as it compresses in the point mass potential estab¬ 
lished by the cold gas. This type of central heating, 
with temperatures extending to 2 — 3 x 10 7 K, has 
not been observed. We have examined the convec¬ 
tive stability of this core of hot gas and find that it is 
stable apart from a small rather insignificant region 
< 500 pc from the very center. However, the massive, 
compact cloud of cold gas that creates the central po¬ 
tential would also increase the central stellar velocity 
dispersion in a manner similar to that of a central 
massive black hole. But the total mass of cooled gas 
in this model, M cg (t n ) = 4.8 x 10 10 M 0 , is about 15 
times more massive than the largest known central 
black hole in elliptical galaxies, 3 x 10 9 M 0 for M87 
(Kormendy & Richstone 1995). We conclude that all 
the interstellar gas that cools since t* cannot reside 
within 2 or 3 kpc of the galactic center since its grav¬ 
ity would unrealistically heat both the stars and hot 
gas. We therefore propose that this cooled gas goes 
into a stellar subsystem or disk of dimension > 3 kpc 
where its gravitational influence is greatly lessened. 
This assumption is consistent with our neglect of the 
gravity of the cold gas and, since M cg (t n ) -C M*, the 
additional contribution of the cold gas mass to stars 
is relatively insignificant. The stellar mass within 
1 kpc M, (1 kpc) = 4.1 x 10 10 M & is comparable to 
M cg (t n ), but the stellar mass within 3 kpc is consid¬ 


erably larger, 13 x 10 10 Mq. Only about half of the 
cold gas mass was available to form additional stars 
several Gyrs after f*, more recently cooled gas must 
be accommodated in another way if there is to be 
no optical evidence for recent star formation. Disks 
are attractive, but we recognize that the formation 
of disks in rotating ellipticals (Brighenti & Mathews 
1996; 1997a) is accompanied by centrally flattened X- 
ray images that have not yet been observed. Options 
for the final disposition of cold gas will be clarified in 
the future when we include galactic rotation in models 
such as those described here. 

To illustrate the critical importance of additional 
circumgalactic gas in determining the global proper¬ 
ties of hot gas in ellipticals, we show in Figure 8 the 
density and temperature variations at time t n for the 
standard model when the dark halo within 150 kpc 
is non-varying and the only source of interstellar gas 
is stellar mass loss. The dramatic inadequacy of this 
model is seen in the density plot where the mass of 
hot gas is far too low in the outer galaxy. But the gas 
temperature is also too low. The central gas temper¬ 
ature in this model is consistent with observed tem¬ 
peratures only at the very center where the observed 
temperature is nearly equal to the stellar virial tem¬ 
perature (T*) « 10 7 K. 

It is important to emphasize that the low temper¬ 
ature in this model cannot be brought into agreement 
with the observations by increasing the current and 
past rate of Type la (or Type II) supernovae. There 
are two reasons for this. First, increasing the SNIa 
rate also increases the interstellar iron abundance far 
above observed values (e.g. Loewenstein & Mathews 
1991). Second, Brighenti & Mathews (1998) found 
that the gas temperature profile at t n is surprisingly 
insensitive to increases in the SNIa rate until that 
rate becomes large enough to drive a galactic wind. 
When a wind develops the gas temperature rises far 
too much and the gas density falls very far below ob¬ 
served values at every galactic radius. We could find 
no value of SNu(f) that corrected the obvious temper¬ 
ature discrepancies like that shown in Figure 8. Davis 
and White (1996) found that gas temperatures in all 
bright ellipticals exceed the stellar virial temperature 
T* by factors of ~1.5, as we find in Figure 3. We 
conclude that the long term accumulation of cosmic 
gas in the dark halo - that naturally results in both 
densities and temperatures similar to those observed 
- is essential to understand the nature of hot gas in 
ellipticals today. 



8. VARIATIONS ON THE STANDARD MODEL 

We now briefly discuss additional models in which 
some of the parameters are varied from those used 
in the standard model; unless otherwise indicated the 
parameters are the same as those used in the standard 
model. 


8.1. Variation of 0/, 

Within the restrictions imposed by primordial nu¬ 
cleosynthesis, flb(H/50) 2 = 0.05 ±0.01 is tightly con¬ 
strained. Nevertheless, we find that our models are 
sensitive to these small changes in Ob. For example 
Figure 9 shows the density and temperature profiles 
at t n = 13 Gyrs for f lb = 0.04 and 0.06. The overall 
agreement with the observed density and tempera¬ 
ture profiles is noticeably degraded compared to the 
fib = 0.05 standard model described earlier. Evi¬ 
dently the amount of cosmically inflowing gas that 
enters the galaxy has a profound influence on its X- 
ray appearance. 

The inflow of gas at the outer edge of the galaxy, 
some of which may have previously flowed out of the 
galaxy during the release of SNII energy, also strongly 
influences the metallicity of hot gas deep inside the 
galaxy. Interstellar gas within the stellar galaxy is 
a blend of gas that has come from stellar mass loss 
and from cosmic inflow; its net metallicity is also a 
blend. To illustrate the contribution from these two 
sources of gas and its dependence on f lb, we show in 
Figure 10 the fraction of gas at time t n that has come 
from stellar mass loss within the galaxy since time f* 
when the SNII energy was released. In the fib = 0.04 
solution almost all of the gas within the galaxy must 
have about the same average metallicity as the stars. 
In fact some of the gas lost from the stars in this 
model since t* has moved out beyond the stellar sys¬ 
tem which extends only to 100 kpc. By contrast, only 
0.7 and 0.3 of the gas at r e = 8.57 kpc has come from 
post-f* stellar mass loss for fib = 0.05 or 0.06 respec¬ 
tively. The blending of two sources of interstellar gas, 
having different metallicities, complicates further the 
current controversy regarding the apparent discrep¬ 
ancy of the gaseous iron abundance with that in the 
stars, ZFe,gas < ^Fe,* (Arimoto et al. 1997; Renzini 
1998). Our fib = 0.04 model may be inconsistent with 
this inequality. 


8.2. Variation of SNII Energy and f* 

The discrepancies seen in Figure 9 must not be 
taken too seriously since some of the misalignment 
with observations can be removed by adjusting the 
total SNII energy (rju) released and the epoch i* of 
galaxy formation. For example, when fib = 0.06 the 
excess gas density in Figure 9 can be expelled from 
the galaxy by increasing the total SNII energy re¬ 
leased at time by 2. Alternatively, if the SNII 
energy is left unchanged, the fib = 0.06 model can 
be improved by forming the stars at an earlier time, 
f* ~ 1.5 Gyr, since there is just enough baryonic gas 
available at that early time to construct the stellar 
galaxy of mass M*. Conversely, the fib = 0.04 solu¬ 
tion in Figure 9 can be improved by lowering the SNII 
energy by 2 and by assuming t* « 3 Gyrs. Several 
models adjusted in this manner are shown in Figure 
11. As t* is increased the mass of gas within the ac¬ 
cretion shock is greater so the post-SNII temperature 
is lower; this follows since we suppose that the SNII 
energy is deposited within the shock radius at f*. In 
all our models we have assumed complete efficiency 
in converting SNII energy to the hot gas, e sn = 1, so 
we have probably somewhat overestimated the SNII 
heating. Each model depends on the product e sn r]u 
so if e sn < 1 then ?;// would need to be increased by 
l/e sn for identical results. 

In principle the degeneracy in the choice of the 
parameters fib, r/n and f* can be partially removed 
by considering in detail the enrichment by elements 
produced in SNII events. We shall not attempt this 
here, but simply note that to compensate a relatively 
small, ±0.01 change in Ob requires a factor of ~ 2 
change in the SNII energy or in the mass available 
to be heated. It is therefore unlikely that successful 
models are possible if Ob is much greater or less than 
the bounds set by nucleosynthesis. 

Finally, we have computed models in which SNII 
heating is ignored altogether. In such models the gas 
density at t n greatly exceeds (by factors of 2 - 5) 
the observed gas density at all galactic radii. The 
amount of cold gas at time t n is also unacceptably 
large, M cg (t n ) > 10 11 M 0 . We conclude that SNII 
heating plays a critical role in establishing the amount 
of hot gas observed today in massive elliptical galax¬ 
ies. 
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8.3. Variation of Cooling or Heating 

In our standard model we allow for a modest 
negative gradient in the gas metallicity by increas¬ 
ing A (T, z) toward the center of the galaxy. When 
compared to models with spatially constant A, the 
dA/dr < 0 assumption slightly improves the fit to 
the observed gas density for r ^ 10 kpc. There 
is an inconsistency with this approach since the ob¬ 
served densities have been determined assuming uni¬ 
form metallicity. However, the sensitivity of our mod¬ 
els to reasonable variations of A is not large. 

We have also explored the hypothesis advanced by 
Tucker & David (1997), Binney & Tabor (1995), and 
Ciotti & Ostriker (1997) that an additional heating 
process is present in the hot gas that reduces the 
total amount of cold gas that forms over a Hubble 
time. However, we have not succeeded in finding a 
heating source term for the thermal energy equation 
that significantly improved our models. If a small 
amount of heating due to some spatially distributed 
source is present, it is either rather ineffective or, 
when increased, generates outward propagating waves 
or shocks that result in peculiar density and tempera¬ 
ture profiles and unsatisfactory X-ray surface bright¬ 
ness distributions. The limiting version of such mod¬ 
els are those in which there is no radiative cooling at 
all, as if heating and cooling were in perfect balance, 
as some have suggested. When the standard model 
is recalculated with A = 0, for example, we find that 
the gas density in r < 40 kpc rises as n oc r~ 2 ' 5 , 
completely diverging from the observations; the gas 
temperature is too low by ~ 0.5 in this same region. 
In general therefore, we find that additional sources 
of distributed heating can introduce undesirable side 
effects that tend to degrade the agreement with ob¬ 
servations. 

8.4. Two kinds of “dropout” 

The deviation between our theoretical models and 
the observed gas density is almost always in the sense 
that the model profile is too centrally peaked. This 
type of deviation has been noticed for many years 
and has led to the hypothesis of distributed mass 
“dropout” in which hot gas is assumed to disappear 
from the flow according to some radial dropout func¬ 
tion that can be varied to obtain agreement with ob¬ 
servations (e.g. Fabian & Nulsen 1977; Sarazin & 
Ashe 1989). Originally it was hoped that the mass 
dropout could be understood as some type of ther¬ 


mal instability resulting from an inhomogeneous flow. 
Unfortunately, detailed hydrodynamical calculations 
have shown that thermal instabilities are disrupted 
by a variety of gas dynamical instabilities before an 
appreciable amount of cold gas forms; see Brighenti 
& Mathews (1998) for a more detailed discussion of 
dropout and additional references. 

Nevertheless, in view of the popularity of mass 
dropout models we illustrate in Figure 12 the den¬ 
sity and temperature structure at time t n in the stan¬ 
dard model when an additional sink term —q(p/tdo) 
is included on the right side of the continuity equa¬ 
tion. The model in Figure 12 is based on q = 1 and 
the dropout time is given by the local radiative cool¬ 
ing time tdo = 5m p kT/2ppA (Sarazin & Ashe 1989). 
The dense, cooling regions that are dropping out con¬ 
tribute substantially to the total emission, reducing 
the apparent gas temperature T 0 b s below that of the 
the undisturbed ambient gas T; both temperatures 
are shown in Figure 12. 

However, almost the same density profile can be 
found if the stellar mass ejection rate a* is uniformly 
lowered. For comparison, a model based on lowering 
a* by 0.5 at all times is shown with dotted lines in 
Figure 12. 

A decrease in a*(t) can be understood from two 
astrophysical arguments. First it is possible to re¬ 
duce a* by altering the (power law) IMF parameters 
x, me and m u . Unfortunately the stellar mass to light 
ratio also changes. For example, for values of x, me 
and m u that reduce a*(t n ) by 2, the mass to light ra¬ 
tio M/Ly must increase by ~ 2 for any choice of the 
three adjustable IMF parameters. Such a large stellar 
M/Ly is inconsistent with the value determined by 
van der Marel (1991) and with the excellent fit to the 
X-ray data for NGC 4472 using his M/Ly (Brighenti 
& Mathews 1997b). Second, a* can be “virtually” 
lowered if some of the gas ejected from evolving stars 
(having a normal IMF) never enters the hot gas phase. 
This idea was discussed by Mathews (1990) but dis¬ 
carded because of the many hydrodynamic instabili¬ 
ties that can disrupt stellar ejecta, increase its surface 
area and promote dissipative and conductive heating 
by the hot gas environment. However, if these diffi¬ 
culties can be avoided, the results in Figure 12 sug¬ 
gests that observations of NGC 4472 may be consis¬ 
tent with this second type of “dropout”. 
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9. TIDAL TRUNCATION AND THE RANGE 
OF X-RAY PROPERTIES AMONG EL¬ 
LIPTICALS 

For many years X-ray astronomers have noted that 
the X-ray luminosities of ellipticals L x span a huge 
range for galaxies of similar optical luminosity Lb 
(Eskridge, Fabbiano & Kim 1995). In spite of much 
effort, no intrinsic property of elliptical galaxies had 
been found that correlated with this scatter (e.g. 
White & Sarazin 1991). However, in a recent paper 
(Mathews & Brighenti 1998) we showed that L x /Lb 
correlates significantly with the physical size of the 
X-ray source, L X /L B oc (r ex /r e )°- 60±0 - 30 , where r ex 
is the radius that contains half of the X-ray lumi¬ 
nosity in projection. The discovery of this correlation 
was not accidental, but was inspired by the hypothesis 
that halo gas and dark matter could be tidally traded 
among massive early type galaxies in small groups. 
Such a tidal exchange of halo material would result in 
large, tidally dominant and centrally located ellipti¬ 
cals having huge, extensive halos and secondary donor 
(giant) ellipticals with more modest, tidally truncated 
dark matter and hot gas halos. Because of the short 
dynamical time scales in galaxy groups, many authors 
have noted that groups are ideal for both the produc¬ 
tion of ellipticals by mergers and for promoting tidal 
warfare among the component galaxies (Merritt 1985; 
Bode et al. 1994; Garciagomez et al. 1996; Athanas- 
soula et al. 1997; Dubinski 1997). It was this hy¬ 
pothesis that led us to seek the correlation between 
L x /L b and r ex /r e . 

We can now test this hypothesis by truncating the 
dark matter and hot gas halos in the models described 
in this paper. Although we cannot describe fully real¬ 
istic tidal effects with our spherically symmetric mod¬ 
els, tidal truncations can be approximately simulated 
simply by removing all dark matter and hot gas be¬ 
yond a truncation radius rtr at some time tt r in the 
past. For this purpose we begin not with a model that 
describes NGC 4472 at time t n , but with an elliptical 
with somewhat higher L x and containing more hot 
gas. This can be accomplished simply by lowering 
the SNII energy below that of our standard model. 
Immediately after the truncation, hot gas just within 
rtr flows outward and a rarefaction wave passes into 
the interstellar gas at the sound speed. After a sound 
crossing time from ry r to the center of galaxy, ~ 1 
Gyr, a new cooling flow equilibrium is established. 

In Figure 13 we show as dashed lines the final den¬ 


sity and temperature profiles at time t n for the stan¬ 
dard model but modified with pu = 0.3 p s td- As ex¬ 
pected, this model is overdense relative to NGC 4472 
at large radii. Also shown with solid lines are the cur¬ 
rent profiles of a model that was truncated at t tr = 9 
Gyrs at radius r tr = 400 kpc but is otherwise identical 
to the first model. The outer density profile of NGC 
4472 is nicely fit by the truncated model and agree¬ 
ment with observed gas temperatures is satisfactory. 
Since the size of the observed X-ray image of NGC 
4472, r ex = 2.87 r e , is average for X-ray bright ellip¬ 
ticals of similar Lb, it is very likely that NGC 4472 
suffered such a truncation in its past. This truncation 
could have been accomplished by nearby group galax¬ 
ies which may not be as massive as 4472. We find that 
the final n(r) and T(r) profiles are rather insensitive 
to the exact time when the truncation occurred. 

Finally, in Figure 14 we show the trajectory in the 
(L x ILb, r ex /r e )-\Amie caused by truncations of this 
pii = 0.3 p s td model at five different radii rt r , all at 
time ttr = 9 Gyrs. Data for the observed galaxies is 
provided in Table 1 of Mathews & Brighenti (1998). 
The large x at the right in Figure 14 shows the final 
location of the untruncated, gas-rich galaxy shown 
in Figure 13. The position of each small x on the 
left shows the final locus of the truncated models. 
Note that models truncated at r tr = 300 and (partic¬ 
ularly) 400 kpc lie rather close to the observed locus 
of NGC 4472, as one would expect from Figure 13. 
As models of gas-rich ellipticals are truncated at dif¬ 
ferent radii, they move precisely along the correlation 
among observed ellipticals. This most encouraging re¬ 
sult strongly supports our hypothesis that the spread 
in the ( L x /Lb , r ea; /r e )-plane is generated by tidal 
exchanges of halo material in small groups. 

10. FURTHER REMARKS 

We have not considered here the detailed metal en¬ 
richment of the hot gas expected from Type la and 
II supernovae; this will be the subject of a future 
study. Clearly, accurate observations and theoreti¬ 
cal predictions of the radial variation of the abun¬ 
dance of iron, silicon and other elements in the hot 
gas will provide important additional constraints on 
the formation and evolution of elliptical galaxies. In 
the absence of this detailed information, it is inter¬ 
esting to compare the total amount of iron produced 
by SNII in NGC 4472 with the iron present in the 
stellar system today. According to Trager (1997) the 
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ratio of iron to a elements is less than solar in most 
bright ellipticals, suggesting that the mean stellar 
iron abundance in NGC 4472 is about ( ZFe ,*) ~ 0.7. 
The total mass of stellar iron in NGC 4472 is then 
Mp e « (zFe,*}Mt,/lA = 6 x 10 8 M & where 1.4 is the 
ratio of total to hydrogen mass. If ~ 0.14 Mg of iron 
is produced in each SNII event (Gibson, Loewenstein, 
& Mushotzky 1998), then Type II supernovae should 
produce a total iron mass of ~ 7 x 10 8 M 0 , compa¬ 
rable but somewhat larger than the mass of iron ob¬ 
served in stars today. Evidently the iron-producing 
capability of the Salpeter IMF is sufficient to account 
for all the iron in stars. In making this estimate we 
have assumed that the iron enrichment due to Type 
la supernovae goes primarily into the hot gas phase 
rather than into stars; the total amount of iron in the 
hot gas within 150 kpc is rather small, ~ 3 x 10 7 Mg.. 

In general the density profiles in our most suc¬ 
cessful models resemble power laws throughout most 
of the observed galaxy while the observed n(r) sug¬ 
gests a double power law with a flatter slope within 
~ 16 kpc. As a result our models predict a steeper 
gas gradient within the central ~ 16 kpc than is ob¬ 
served. It is this small deviation that the “dropout” 
hypothesis was designed to correct. While future ob¬ 
servations will clarify this curious systematic devia¬ 
tion, it is possible that new physical effects become 
important in r ^ 10 kpc. We have already noted 
that a small, but realistic galactic rotation could flat¬ 
ten the azimuthally-averaged gas density profile in 
this region. Another possibility within ~1 kpc is 
the growth of magnetic pressure toward the galac¬ 
tic center. Brighenti & Mathews (1997b) have shown 
that the total mass predicted from hydrostatic models 
is systematically lower than the known stellar mass 
within the central kiloparsec in each of the three 
bright ellipticals studied with Einstein HRI. This de¬ 
viation can be understood if the gas pressure is as¬ 
sisted by an additional component of non-thermal 
pressure. The required magnetic field is large - we 
estimated B ~ 10~ 4 gauss for NGC 4636 but not 
unexpected. Mathews & Brighenti (1997) showed 
that fields of this magnitude are a natural outcome 
of the ejection of fields by mass-losing stars. Fields 
ejected in stellar envelopes are amplified by turbulent 
dynamo action in the interstellar medium and by the 
natural intensification that accompanies the inward 
flow of interstellar gas in which the magnetic pressure 
increases more rapidly than the gas pressure. The im¬ 
portance of magnetic or other non-thermal pressure 


components in the interstellar cores of ellipticals will 
be clarified when AXAF observations become avail¬ 
able. 

11. CONCLUSIONS 

We have described here the overall evolution of 
massive elliptical galaxies with particular attention 
to the hot interstellar gas. Although the stellar sys¬ 
tems in giant ellipticals are thought to be very old, the 
mixing of the stellar orbits masks many details of the 
dynamical and enrichment processes that occurred 
during the time of galaxy formation. The metallic- 
ity gradient in the stars provides some clue to these 
early events but its interpretation is far from obvious. 

In contrast, the hot interstellar gas in ellipticals to¬ 
day does retain a large amount of information about 
star formation and metal enrichment events that oc¬ 
curred during the earliest phases of galaxy evolution. 
Interstellar gas in ellipticals can help us solve one of 
the thorniest problems in modern studies of galaxy 
formation: the importance of “feedback” SNII energy 
on the gas. Among the major feedback uncertainties 
are: the total amount of supernova energy involved, 
its efficiency in heating the gas, and the time that 
SNII enrichment occurred relative to the dynamical 
formation of the galaxies. Regarding this last point, 
it is currently unclear if some fraction of the metals 
in the hot gas were made and ejected from an earlier 
generation of low mass galaxies before the formation 
of galaxy groups, or if all the metals were made dur¬ 
ing the formation of group galaxies that merged soon 
thereafter to create the ellipticals observed today. 

We believe that studies of the creation and evolu¬ 
tion of hot gas in large ellipticals provide a unique 
insight into the feedback problem and may help us 
unravel some of these complex events that occurred 
so long ago. Many lines of astronomical evidence in¬ 
dicate that the formation sites of most massive el¬ 
lipticals are small groups of galaxies that are very 
old. Dynamical studies of the evolution of galax¬ 
ies in elliptical-forming groups indicate that ellipti¬ 
cals formed early in the dynamical evolution of these 
groups when mergers and tidal interactions were most 
likely; if so, both the elliptical galaxy and the hot in¬ 
terstellar gas surrounding it are very old. Remark¬ 
ably, the hot gas we see today still contains a record 
of these past events. Since SNII enrichment is accom¬ 
panied by an increase in entropy due to heating of the 
gas, enrichment and gas dynamics are closely linked. 
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In this study we have described the first compre¬ 
hensive theoretical description of the evolution of hot 
gas in massive ellipticals that is fully consistent with 
important cosmological constraints. The goal of our 
calculations is not merely to explain the total mass 
and energy in the interstellar medium, but to account 
for the detailed radial dependence of these quantities 
throughout the observable range. In spite of the sim¬ 
plicity of our galaxy-forming model and our choice 
of a currently unfashionable flat cosmology, we have 
been remarkably successful in explaining X-ray obser¬ 
vations of bright ellipticals in a cosmological context 
that conserves baryonic and dark mass, allows for all 
energy sources involved, and includes the detailed in¬ 
teraction of gas with mass-losing stars. 

We now summarize some of the main conclusions 
of the calculations we have discussed here: 

(1) Hot gas density and temperature profiles in 
massive ellipticals cannot be understood solely with 
gas ejected from the old stellar population that ac¬ 
counts for the optical light observed today. 

(2) Low interstellar gas temperatures in models 
without cosmic inflow cannot be heated to the ob¬ 
served values by current and past Type la supernovae; 
if this is attempted, strong galactic winds develop 
that are incompatible with observation. 

(3) Excellent agreement with observed hot gas 
mass and thermal energy distributions - n(r) and 
T(r) - can be achieved by the addition of cosmic 
gas from secondary infall. This gas shocks and com¬ 
presses to the virial temperature of the galactic halo, 
accounting for the hot gas temperatures observed in 
all massive ellipticals. As this gas flows deeper into 
the galaxy, it mixes with somewhat cooler gas ex¬ 
pelled from the stars, resulting in positive gas temper¬ 
ature gradients ( dT/dr > 0) in r ^ 3r e as observed in 
most or all bright ellipticals (see Fig. 1 of Brighenti 
& Mathews 1997b). 

(4) Within the framework of a flat cosmology, we 
can account for observed gas density and temperature 
profiles only for a rather narrow range in baryon mass, 
H;, « 0.05 ± 0.01, the same range that is consistent 
with primordial nucleosynthesis. 

(5) The fraction of gas near the optical effective 
radius that derives from stellar mass loss is strongly 
and inversely dependent on fib when other parameters 
are held constant. 

(6) However, our successful fits to observed n(r) 
and T(r) are degenerate in some of the parameters. 


For example, increasing fit, is approximately equiv¬ 
alent to (i) decreasing the SNII energy (??//) or its 
efficiency in heating the gas (e sn ) or (ii) increasing 
the time f* that the bulk of galactic stars formed so 
that more gas (within the accretion shock) needs to 
be heated by SNII. Both (i) and (ii) reduce the en¬ 
ergy per gram delivered to the gas by SNII explosions. 
This degeneracy is due in part to the simple, single 
burst assumption we have made in forming stars and 
our assumption that SNII energy is deposited within 
the shock radius at time f*. More complex mod¬ 
els that consider the details of metal enrichment and 
other options for the spatial distribution of SNII en¬ 
ergy may remove some or all of this degeneracy. 

(7) The current fraction of mass in baryons is ex¬ 
pected to have a large minimum surrounding large el¬ 
lipticals and their associated groups. This minimum 
results from an outward displacement of baryons that 
occurred when the energy of Type II supernovae was 
released. As a result, observational measurements of 
flb/f1 within large galactic radii can either overesti¬ 
mate or underestimate the true cosmic value of this 
ratio. 

(8) Since the epoch of star formation the amount 
of gas that has cooled is comparable to the total 
mass of hot gas observed today. This gas cannot 
all cool within the central few kiloparsecs since its 
mass would unrealistically heat both the gas and stars 
there. However, even a tiny galactic rotation is suf¬ 
ficient to avoid this central concentration of cooled 
gas. 

(9) Finally, when a relatively gas-rich version of 
one of our models is ticlally truncated at some time 
in the past, the X-ray luminosity L x and the effective 
radius for X-ray emission r ex both decrease. The vec¬ 
tor due to truncation in the ( L x /Lb , r ea ,/r e )-plane 
lies almost exactly along the correlation discovered 
by Mathews & Brighenti (1998). This give addi¬ 
tional support to the notion, described by Mathews 
& Brighenti (1998), that the enormous spread in L x 
for fixed galactic Lb is due in some large measure to 
the tidal exchange of halo gas and dark matter among 
galaxies within groups of galaxies. 
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light ratios for us. Our work on the evolution of hot 
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Fig. 1.— Upper panel: Open circles show ROSAT ob¬ 
servations of plasma temperature (upper curve) and 
density (lower curve). Filled circles are derived from 
Einstein HRI observations. Lower panel: Solid curve 
is the total gravitating mass based on assuming hy¬ 
drostatic equilibrium; Dotted line is a de Vaucouleurs 
fit to the stellar core; Long dashed line is an NFW 
dark halo that fits the outer galaxy. Dot-dashed line 
is the sum of the de Vaucouleurs and NFW mass pro¬ 
files. 

Fig. 2.— NFW dark halo (solid line) and mass in the 
secondary inflow (dashed lines) shown at four times: 
0.5, 3.0, 8.0 and 13.0 Gyrs from left to right. At large 
radii the mass increases as r 3 . 

Fig. 3.— Upper panel: gas density profile of the stan¬ 
dard model compared to observations of NGC 4472; 
Lower panel: gas temperature profile of the standard 
model compared to observations of NGC 4472. 

Fig. 4.— Flow variables in the standard model at 
three times: 4 Gyrs (dotted lines), 9 Gyrs (dashed 
lines) and 13 Gyrs (solid lines). Upper left: gas den¬ 
sity, Upper right: gas pressure, Lower left: gas veloc¬ 
ity, and Lower right: gas temperature. 

Fig. 5.— Solid line: Variation of the baryon fraction 
at time t n = 13 Gyrs in the standard model. Dashed 
line: The baryonic contribution of the gas. As the 
radius increases the fraction becomes f If, = 0.05. 

Fig. 6.— History of various mass components in 
the standard model. The stellar system is created at 
time i* = 2 Gyrs. Solid line: Variation of the mass 
of hot gas within the accretion shock radius, Dashed 
line: Variation of the mass of hot gas within 150 kpc, 
Dotted line: Variation of cold gas (T < 10 5 K), Dot- 
dashed line: Variation of the mass in stars. 

Fig. 7. Density and temperature profiles for the 
standard model (as in Fig. 3) but including the grav¬ 
ity of the cold gas, assumed to be a central point mass. 
Upper panel: gas density profile compared to obser¬ 
vations of NGC 4472; Lower panel: gas temperature 
profile compared to observations of NGC 4472. 

Fig. 8.- Density and temperature profiles for the 


standard model (as in Fig. 3) but with stellar mass 
loss as the only source of interstellar gas. Upper panel: 
gas density profile compared to observations of NGC 
4472; Lower panel: gas temperature profile compared 
to observations of NGC 4472. 

Fig. 9. Density and temperature profiles for models 
similar to he standard model in Figure 3 but with 
Lib = 0.06 (solid line) and fit, = 0.04 (dashed line). 

Fig. 10.— Fraction of gas ejected from galactic stars 
since time i* for the standard model with Lf, = 0.05 
(solid line), and for similar models computed with 
Lib = 0.04 (dashed line) and Lib = 0.06 (dotted line). 

Fig. 11. Several models with various values of Lib, 
f* and t)ii . Solid line: Lib = 0.05, t* = 2 Gyrs and 
Vn = O.Grjstd- Dashed line: Lib = 0.04, t,„ = 3 Gyrs 
and 7ju = O.bijstd- Dotted line: Lib = 0.06, f* = 2 
Gyrs and p n = 2 r) st d- 

Fig. 12.— Density (upper panel) and temperature 
(lower panel) distributions of dropout and low a* 
models compared with observations of NGC 4472. 
Dashed line: Dropout in the standard model with 
q = 1, Solid line: Dropout in standard model with 
reduced SNII energy, rju = 0.6rj st d, In the dropout 
models the upper (thin) lines show the temperature 
variation of the background flow and the lower (thick) 
line show the apparent temperature T 0 b s allowing for 
the cooling regions. Dotted line: Standard model with 
a* reduced by a factor 0.5 at all times and galactic 
radii. 

Fig. 13.— Dashed line: Density and temperature 
profiles at t n for a model similar to he standard model 
in Figure 3 but with pu = 0.3 r/ s td- Solid line: A 
repeat calculation of the same model at t n which was 
truncated at radius rt r = 400 kpc at time tt r = 9 
Gyrs. 

Fig. 14.— A plot of the (L x /Lb, r ea; /r e )-plane show¬ 
ing the correlation among observed ellipticals. The 
large x at the right is the locus at t n of the untrun¬ 
cated model shown in Figure 13. The small xs that 
extend away from the untruncated solution toward 
the lower left show the locus at t n of models trun¬ 
cated at time tt r = 9 Gyrs at five decreasing trun- 
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cation radii r tr = 500, 400, 300, 200, and 100 kpc. 
As rtr decreases the locus of the final position in the 
(. L x /Lb , r e x/r-e)-plane moves from the untruncated 
solution almost exactly through the observed correla¬ 
tion. 
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